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I. INTRODUCTION 

The influence of impurities, defects, or other types 
of quenched disorder on the properties of phase tran- 
sitions has aroused the interest of physicists for more 
than three decades (see Ref. 1 for an overview of some 
of the early work). Recently, this field has reattracted 
considerable attention as it has become clear that dis- 
order effects are generically much stronger at zero- 
temperature quantum phase transitions than at classi- 
cal thermal phase transitions. This leads to unconven- 
tional phenomena such as power-law quantum Griffiths 
singularities, 2-4 infinite-randomness critical points with 
exponential rather than power-law scaling, 5,6 or even 
smeared phase transitions. 7,8 A recent review of part of 
this physics can be found in Ref. 9. 

The main reason for the enhanced disorder effects at 
quantum phase transitions is that the disorder is per- 
fectly correlated in imaginary time direction. Because 
imaginary time acts as an extra dimension at a quantum 
phase transition (and becomes infinitely extended at zero 
temperature), one is effectively dealing with defects that 
are "infinitely large" in this extra dimension. Thus, they 
are much harder to average out than conventional finite- 
size defects. 

This implies that similarly strong effects can be ex- 
pected at a classical thermal phase transition if the dis- 
order is perfectly correlated in one or more space dimen- 
sions. Indeed, it has been known for a long-time that 
the McCoy- Wu model, a classical two-dimensional Ising 
model with disorder perfectly correlated in one of the 
two dimensions, exhibits an unusual phase transition. 
In a series of papers, 10-13 McCoy and Wu developed a 
transfer-matrix approach to this model and showed that 
the specific heat is smooth across the ferromagnetic phase 
transition while the susceptibility is infinite over an en- 
tire temperature range. Fisher 5,6 later achieved an es- 
sentially complete understanding of this transition by 



means of a strong-disorder renormalization group (us- 
ing the equivalence between the McCoy- Wu model and 
the one-dimensional random transverse-field Ising chain) . 
He found that the critical point is of infinite-randomness 
type, and it is accompanied by strong power-law Griffiths 
singularities. Largely due to the fact that the McCoy- Wu 
model is difficult to realize in nature, these predictions 
have (to the best of our knowledge) not been experimen- 
tally verified yet. 

In this paper, we present another classical system ex- 
hibiting an exotic infinite-randomness critical point, viz., 
a randomly layered three-dimensional Heisenberg mag- 
net. This system is more easily realizable in experiment 
than the McCoy- Wu model as it can be produced by de- 
positing random layers of two different ferromagnetic ma- 
terials. Moreover, because of its three-dimensional char- 
acter, it permits bulk thermodynamic measurements. We 
investigate the phase transition in this model by means 
of a strong-disorder renormalization group which allows 
us to determine the critical behavior exactly. 

Our paper is organized as follows: In Sec. II, we in- 
troduce the randomly layered Heisenberg model and give 
heuristic arguments for the strong disorder effects. In 
Sec. Ill we explain our theoretical approach. Results on 
the thermodynamics are given in Sec. IV, while the ex- 
perimentally important finite-size effects are discussed in 
Sec. V. Section VI is devoted to the dynamical behavior 
at the phase transition. We conclude in Sec. VII. 



II. RANDOMLY LAYERED HEISENBERG 
MODEL 

We consider a three-dimensional Heisenberg ferromag- 
net consisting of a random sequence of layers made of 
two different ferromagnetic materials, as sketched in Fig. 
1. This system can be modelled by a classical Heisenberg 
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FIG. 1: Schematic of the layered magnet: layers of two dif- 
ferent ferromagnetic materials are arranged in a random se- 
quence. 

SO WO WD SD 
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FIG. 2: Schematic phase diagram of the randomly layered 
Heisenberg magnet (1). SD and SO denote the conventional 
strongly disordered and strongly ordered phases, respectively. 
WD and WO are the weakly disordered and ordered Griffiths 
phases. T c is the critical temperature while T u and 21 mark 
the boundaries of the Griffiths phase. 

Hamiltonian on a cubic lattice given by 
H = — jj (S r -S r+ x+S r -S r+ y)— J^~ S r -S r+ £. (1) 

r r 

Here, S r is a three-component unit vector on lattice site 
r, and x, y, and z are the unit vectors in the coordinate 
directions. The exchange interactions within the layers, 
jj, and between the layers, J^~, are both positive and 
independent random functions of the perpendicular co- 
ordinate z. 

To develop a heuristic understanding of the randomly 
layered Heisenberg model, we first consider the case of 
all being identical, J^r = J x , while the jj are drawn 
from a binary probability distribution 

P(jH) = (1 -p)5(jH - J u ) + pS(J - J) (2) 

with J u > J; . Here, p is the concentration of the "weak" 
layers while 1 — p is the concentration of the "strong" 
layers. More general distributions will be considered in 
the next section. 

Let us now discuss the behavior of the model (1) qual- 
itatively (see Fig. 2). At sufficiently high temperatures, 
the system will be in a conventional (strongly disordered) 
paramagnetic phase with a finite magnetic susceptibility 
which increases upon lowering the temperature. Below 
a temperature T u (which is the transition temperature 
of a hypothetical system containing strong layers only, 



Jz = J u ), rare thick slabs of strong layers develop local 
order while the bulk system is still nonmagnetic. This 
is the weakly disordered Griffiths phase. The Griffiths 
phase continues below the actual critical temperature 
T c down to a temperature Ti (which is the transition 
temperature of a hypothetical system containing weak 
layers only, jj = J;). In the weakly ordered Griffiths 
phase, bulk magnetism coexists with locally nonmagnetic 
slabs. Finally, below T/, the system is in a conventional 
(strongly ordered) ferromagnetic phase. 

To estimate the strength of the Griffiths singularities in 
this system, we need to compare the probability of finding 
a thick slab of strong layers with the contribution such a 
slab can make to thermodynamic quantities such as the 
susceptibility. Simple combinatorics gives the probability 
for finding a slab of L RR consecutive strong layers to be 

w(L RR )~(l-p) L ™ =e-P L ™ (3) 

with p = — ln(l — p). Each such slab is equivalent to a 
two-dimensional Heisenberg model with an effective ex- 
change interaction L RR J U . Because the two-dimensional 
Heisenberg model is exactly at its lower critical dimen- 
sion, the susceptibility of the slab increases exponentially 
with the effective interaction, 9,14 

x(Lrr) ~ e bL ™ (4) 

where b increases with decreasing temperature. The same 
result also follows from a renormalization group analysis 
of the corresponding nonlinear sigma model at its low- 
temperature fixed point 15 or from an explicit large-iV 
calculation (as shown in the next section). 

Thus, the exponential decrease of the rare region prob- 
ability w(L RR ) with size L RR is compensated by an ex- 
ponential increase in the contribution it makes to the 
susceptibility. The total rare region susceptibility in the 
weakly disordered Griffiths phase is obtained by simply 
summing over the contributions of the individual rare 
regions. Up to pre-exponential factors, this yields 

XRR = J dL RR e^ L - . (5) 

The total rare region susceptibility thus diverges once 
b becomes larger than p. Other observables can be 
discussed along the same lines. Equations (3), (4), 
and (5) are analogous to the corresponding relations for 
the McCoy- Wu model 10-13 or (after quantum-to-classical 
mapping) to those of the random transverse-field Ising 
model. 2-4 This suggests that the phase transition in our 
model displays unconventional behavior. In the next sec- 
tion we investigate this question in detail by means of a 
renormalization group method. 

III. STRONG-DISORDER RENORMALIZATION 
GROUP 

In this section we study the ferromagnetic phase transi- 
tion of the randomly layered Heisenberg model by means 
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of a strong-disorder renormalization group. 16,17 Our im- 
plementation of this method follows a recent study of 
dissipative quantum phase transitions. 18,19 We therefore 
only outline the major steps of the calculation, details 
can be found in Ref. 19. 



A. Order parameter field theory 

Our starting point is a Landau-Ginzburg- Wilson 
(LGW) order parameter field theory for an iV-component 
order parameter </>(r). In the absence of disorder, the free 
energy functional reads 



S = j d?r [So 4> 2 {v) + ll{dA{v) f + u<p 4 (r)] . 



(6) 



Here, 6q is the bare distance from criticality, 70 is the 
bare length scale, and u is the </> 4 coefficient. In the pres- 
ence of our layered disorder, 60, 70, and u become random 
functions of the ^-coordinate (the coordinate perpendic- 
ular to the layers) and the derivative term will generally 
be anisotropic. In order to apply the real-space-based 
strong-disorder renormalization group, we discretize the 
continuum LGW theory (6) in the z-dircction, but not in 
the xy-plane. For simplicity, we first consider the large- N 
limit of our LGW theory which allows us to perform all of 
the following calculations explicitly. We will later show 
that the resulting critical point is the same for all N > 2 
including the physically relevant Heisenberg case N = 3. 
The discrete large- N order parameter field theory reads 

s = 5>, + K +7»V)l^(q)l 2 - E ^<Mq)<W-q) 



(7) 

where q is a two-component vector describing the xy- 
momentum. The Lagrange multipliers X z enforce the 
large- N constraints {{(j)^) 2 ) = 1 for the fc-th order pa- 
rameter component in layer z; they have to be deter- 
mined self-consistcntly. The renormalized local distance 
from criticality in layer z is given by e z = S z + X z . In 
the disordered phase, all e z > 0. For the case of a single 
layer, the LGW theory (7) can be solved immediately, 
giving 



7fA 2 e' 



-47T7f /a 2 



(8) 



with A being a momentum cutoff and a the lattice con- 
stant. 



B. Recursion relations 

The basic idea of the strong-disorder renormalization 
group is to successively integrate out local high-energy 
degrees of freedom. In the LGW theory (7), the compet- 
ing local couplings are the local distances from criticality 
e z and the interactions J z . In the bare theory, they are 
independent random variables with distributions Ro{e) 



and Po(^ J ")j respectively. The method relies on these 
distributions being broad and becomes exact in the limit 
of infinitely broad distributions. We will verify this con- 
dition a-posteriori. 

In each renormalization group step, we choose the 
largest local coupling f2 = max{e z , J^}. If it is a distance 
from criticality, say £2 , the unperturbed part of the free 
energy is So = ]C q (e2+72 c L 2 )l ( Mq)| 2 - The coupling of </> 2 
to the neighboring layers, Si — — X) q [^I L< M c l) < / , 2(— q) + 
J^ L 2 (q)?! , 3(— q)], is treated pcrturbatively. Keeping only 
the leading long- wavelength terms that arise in 2nd order 
of the cumulant expansion, we obtain new interactions 

s = -E q ^i L <Ai(q)^3(-q) with 



J 1 - T 1 - 



£2 



(9) 



At the end of the renormalization group step, <j>2 is 
dropped from the action. 

If the largest local energy is an interaction, say , 
we solve the two-layer problem So = SqEz=2 3fe 

7 2 q 2 )l^ z (q)| 2 - E q * / 2 L< Mq)</'3(-q) exactly while treat- 
ing the interactions with the neighboring layers as per- 
turbations. For J '2 >• £2, £3, the two fields and fa 
are essentially parallel; thus they can be replaced by a 
single field cf>2 with an effective renormalized free energy 
functional S = £ q (£~ 2 + 7 2 V ) I fo (q) | 2 . 

After a straightforward but somewhat lengthy 
calculation, 19 the effective distance from criticality of the 
combined layer comes out to be 



£2 = 2 



£2£3 



(10) 



while the length scale parameter renormalizes as 7! = 
72 + 7 2 ■ The new field represents a layer with effective 
moment per site 



M2 = M2 + A*3 • 



(11) 



The interactions of the combined layer with the neigh- 
boring layers are not renormalized. The net result of the 
renormalization group step is the elimination of one layer 
and the reduction of the energy scale ft. 

The structure of the renormalization group recur- 
sion relations (9-11) is identical to those of the one- 
dimensional random transverse-field Ising model 5 ' 6 as 
well as the dissipative quantum rotor model. 18,19 Con- 
sequently (and somewhat surprisingly), the thermal 
phase transition in our randomly-layered classical three- 
dimensional Heisenberg model belongs to the same uni- 
versality class as the quantum phase transitions in the 
one- dimensional random transverse-field Ising model and 
the dissipative quantum rotor chain. 

At first glance, this result seems to suggest that cru- 
cial system characteristics such as order parameter sym- 
metry and dimensionality are rendered unimportant by 
the strong-disorder renormalization group. However, the 
physics turns out to be more subtle. The fact that 
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our randomly layered Heiscnbcrg model and the random 
transverse-field Ising chain are in the same universality 
class is due to a nontrivial interplay between the order 
parameter symmetry and the defect dimensionality. We 
will discuss this point in more detail in Sec. VII in the 
context of a general classification of phase transitions in 
the presence of disorder. 

C. Fixed points 

The renormalization group step outlined in the last 
subsection does not change the lattice topology because 
we remove a full layer in each step. Moreover, the sur- 
viving e and J remain statistically independent. The 
theory can therefore be formulated in terms of individ- 
ual probability distributions P(J _L ) and R(e). Fisher de- 
rived flow equations for these distributions and solved 
them analytically. 5,6 They have three three kinds of non- 
trivial fixed points representing the weakly ordered and 
disordered Griffiths phases as well as the critical point in- 
between. At the critical fixed point, the relative width of 
the distributions P(J _L ) and R(e) diverges, justifying the 
method and giving the critical point its name, infinite- 
randomness critical point. 

The critical behavior is characterized by three expo- 
nents, v = 2, V = 1/2, and <p = (1 + y/5)/2. The expo- 
nent v controls how the perpendicular correlation length 
£_L diverges as the critical point is approached 

u - \sr ■ (12) 

characterizes the spatial correlations perpendicular to 
the layers (in z-direction). 5 is the fully renormalized 
distance from the critical point; it is given by S ~ [ln(e) — 
ln(J ± )]o in terms of the bare variables ([-]o denotes the 
average over the bare disorder distributions.) 

The exponent ip (which is sometimes called the tun- 
neling exponent because of its meaning in the quantum 
problem of Refs. 5,6) relates the perpendicular correla- 
tion length £j_ and the correlation length £|| within the 
layers. The scaling is highly anisotropic, 

M£||/«)~tf ■ (13) 

tp also controls the density nn of layers surviving at en- 
ergy scale f2 in the renormalization procedure. The scal- 
ing form of this variable is given by 19 

nn(6) = [ln(n I /n)}- 1 ^X n [6^ln(n I /n)] , (14) 

where £1/ is a constant of the order of the initial (bare) 
value of fi. The scaling function behaves as A„(0) = 
const and X n (y — > oo) ~ y 1 /^e~ cy where c is a con- 
stant. As a result, the layer density decreases as nn <~ 
[ln(0/ /f2)] -1 /^ at criticality while it behaves as nn <~ 
8 ,J Vl 1 ' z in the disordered Griffiths phase (S > 0). The 
nonuniversal exponent z varies as z <~ in the Grif- 

fiths phase. 



The exponent <p determines how the typical moment 
/zn per site of a surviving layer depends on the energy 
scale f2. The scaling form of fin reads 

^a(5) = [ln(0 7 /n)]^[^ln(0 7 /0)] . (15) 

The scaling function behaves as A^(0) = const and 
A M (y — > oo) <~ y 1 ^^. Thus, at criticality the typical 
moment increases as /in ~ [M^W^)]^ while it behaves 
as nn ~ S^-^ln^j/fi) in the disordered Griffiths 
phase. 

IV. THERMODYNAMICS 

The overall strategy 6 for computing the behavior 
of thermodynamic observables consists in running the 
strong-disorder renormalization group from the initial en- 
ergy scale fti down to the energy scale set by an external 
perturbation such as a magnetic field. The high-energy 
degrees of freedom eliminated in this way do not make 
significant contributions to the long-wavelength physics. 
The surviving layers are very weakly coupled and can be 
treated as independent. In this section we show that the 
resulting thermodynamic behavior of our system is sim- 
ilar to that at an infinite-randomness quantum critical 
point. However, there are significant differences due to 
the fact that we are dealing with a thermal (classical) 
phase transition. 



A. Single-layer results 

We start by considering a single layer with effective 
moment n per site in an external magnetic field (i.e., this 
layer is the result of combining u original layers during 
the renormalization group). The free energy functional 
is given by 

5 ; = £(5 + A + 7 2 q 2 )l^(q)| 2 -nJWN (16) 
q q 

where h(q) is the Fourier transform of the external field 
at wave vector q. This theory is Gaussian, thus the 
partition function and free energy are easily calculated. 
For a uniform magnetic field h, the free energy reads 
Fi(h) = E q Me + 7 2 q 2 ) - £|/i 2 ^ 2 /4e. Here e = 8 + A as 
before, and is the linear size of the layer. The value 
of the Lagrange multiplier A follows from the large-A 
constraint 

2 1 dFx 1 ^ 1_ m 2 ^ 2 , n7 , 

^^q^ = q\7T^ + ^ = l - (17) 

For small fields, fih <C e(h = 0), the first term in the sum 
dominates, yielding e(h) — e(0) + 0(h 2 ) with e(0) given 
by (8). In the opposite limit, [ih ^> e(h = 0), the second 
term dominates, resulting in e(h) = jih/2. 
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The magnetization of the single layer is easily com- 
puted by taking the appropriate derivative of the free 
energy 

mi = -{l/Lf)(8Fi/dh) e = p 2 h/2e ; (18) 

and the zero-field uniform susceptibility is given by 

Xi = /i 2 /2e(0) . (19) 

Other observables can be computed in an analogous fash- 
ion. For instance, the local susceptibility xi,loc takes the 
same form as (19) with p, 2 replaced by p. 

B. Critical point and weakly disordered Griffiths 
phase 

We now combine the single-layer observables with the 
strong-disorder renormalization group results for the den- 
sity (14) and the moment (15) of the surviving layers. In 
the present subsection we focus on the critical point and 
the disordered Griffiths phase while the ordered Griffiths 
phase will be addressed in Sec. IV C. 

The total magnetization in a magnetic field h can be 
obtained by running the renormalization group to the 
energy scale tth = Pn h h. All the surviving layers have 
e <C ph and are thus fully polarized. The total magneti- 
zation per site thus reads 

m(5,h) = nn h (5)nn h (5) 

= Mtoi/tohtf-Wernifi^Mtoi/toh)] (20) 

The scaling function is given by Q m (y) = X n {y)X^{y). 
Now, using the fact that Q,h — Pn h h, we find m <~ 
[ln(fij/ft)]^ _1 /^ (with double- logarithmic corrections) at 
criticality, 5 — 0. This implies that the critical isotherm 
exponent 5 (commonly defined via m <~ h 1 /*) is for- 
mally infinite. In the Griffiths phase, S > 0, we obtain 
m ~ } l V*8 v +»i>Q-m+V*) [ln(0 7 //i)] l+1/z . As long as 
z ~ 5~ v ^ is larger than one (i.e., sufficiently close to 
the critical point) this contribution dominates the regu- 
lar linear-response term. We thus find a non-universal 
power-law singularity in a finite temperature interval 
around the critical point. 

In the zero-field limit, the uniform susceptibility x = 
dm/dh ~ ft, 1 / 2-1 consequently diverges not just at the 
critical point but for all z > 1, again in an entire tem- 
perature range around the critical point. This result can 
also be obtained by summing (19) over all layers using 
the spectral density p(e) = dnn/drt\a= e . In the Griffiths 
phase this gives the following rare region contribution to 
the susceptibility 

x(/w o) ~ { ; z<l ; . (2D 

The specific heat can be obtained by summing the 
single-layer free energy Fi over the spectral density p(e) 
and taking the appropriate derivatives with respect to 
the reduced temperature. As in the McCoy- Wu model, 
the resulting specific heat is smooth across the transition. 



C. Weakly ordered Griffiths phase 

While the order parameter symmetry does not play 
a significant role on the disordered side of the critical 
point where all conventional (non rare-region) excitations 
are gapped, it becomes important on the ordered side of 
the transition where gapless excitations exist even in the 
absence of our rare region physics. This leads to some 
minor differences between our results and those of the 
McCoy- Wu model. 

To determine the spontaneous magnetization on the 
ordered side of the transition, we follow the strong- 
disorder renormalization group flow from 0/ towards 
Q = 0. For small but nonzero \5\, i.e., close to the critical 
point, the flow initially follows the critical trajectory until 
the renormalization group length scale reaches the cor- 
relation length £_l ~ |<5| _ly - This occurs at an energy 
given by ln(fij/%) ~ £f ~ \S\-^ . Beyond this scale, 
the system is essentially ordered, and almost no layers 
will be removed under further action of the renormaliza- 
tion group. 

We can therefore find the spontaneous magnetization 
by counting how many of the original layers survive at 
length scale This leads to 

m ~ na ( IHi ( - Hfij/fi*)]*- 1 /* ~ l^ 1 "^ . (22) 

The order parameter critical exponent thus takes the 
value (3 = v{\ — <fnp). In a small magnetic field h, the 
magnetization picks up a nonanalytic correction which 
can be computed following the methods of Ref. 6. We 
find 

m(h) - m(0) - , (23) 

implying that the (longitudinal) susceptibility x — 
dm/dh - diver ges in the zero-field limit every- 

where in the weakly ordered Griffiths phase. (The trans- 
verse susceptibility is infinite everywhere in the ordered 
phase simply because of the continuous order parameter 
symmetry.) 

Another important property of the ordered phase of a 
continuous symmetry magnet is the spin-wave stiffness 
which can be defined via the change of the free energy 
with a twist in the boundary conditions. In our system, 
we must distinguish the parallel spin- wave stiffness from 
the perpendicular one. To find the parallel spin-wave 
stiffness pi, we apply boundary conditions at x — and 
x = Lh such that the spins at the two ends are at a 
relative angle 9. In the limit of small and large L||, 
the free energy density / depends on 6 as 

/(©)-/(0) = ^l(^) 2 (24) 

which defines pi. 

In our system, the free energy cost due to the twist is 
simply the sum over all layers participating in the long- 
range order. Each layer has the same twisted boundary 
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conditions and the perpendicular bonds (which are not 
twisted) do not contribute. The bare stiffness of a single 
layer is given by 7 2 . Because -f 2 is additive under the 
strong-disorder renormalization group, pi behaves like 
the layer moment per site, pi <~ p. The calculation of 
the total parallel spin-wave stiffness thus proceeds anal- 
ogously to the total magnetization yielding 

p\ ~ 7o 2 = 7o 2 ■ (25) 

If a global twist is applied perpendicular to the layers, 
i.e., between the bottom (z = 0) and top (z = Lj_) of the 
stack, the local twist Z between layers z and z+1 will 
vary from layer to layer according to the local J^r . The 
total free energy cost can be written as 

/(©)-/(0)~^l> e ' (26) 

z 

with p z <~ J z . Minimizing /(©) — /(0) under the con- 
straint J2 Z ®z = ® gives Z <~ 1/ p z and 

/(©)-/(o)~ (^X> _1 ) • ( 2? ) 

To obtain an upper bound for /(O) — /(0), we estimate 
Ez/^ 1 b y its largest contribution, p m *„ - In 
the weakly ordered Griffiths phase, the fixed-point dis- 
tribution of J is gapless, 6 and J^ in vanishes as i]] z 
in the thermodynamic limit L± — > 00. We conclude 
/(O) — /(0) ~ LJ Z_1 , implying that the global perpen- 
dicular stiffness vanishes, p^ = (for z > I). The weakly 
ordered Griffiths phase is thus very peculiar because the 
system displays long-range ferromagnetic order but it has 
no (perpendicular) spin-wave stiffness. 

V. FINITE-SIZE EFFECTS 

The results in Sec. IV were for an infinite system (in 
the thermodynamic limit). Here, we briefly discuss the 
effects of a finite system size in either parallel or perpen- 
dicular direction. 

We start with a finite in-plane (parallel) size L\\. It 
plays the same role as a finite temperature in the quan- 
tum phase transitions in Refs. 5,6,18,19 where the in- 
verse temperature is the system size in imaginary time 
direction. Solving the large-iV constraint for a single 
layer of linear size gives e(£||) = e(oo) + 0(1/ L 2 ) for 
e(oo) 3> V-^m- Here, e(oo) is the thermodynamic limit 
result given in (8). In the opposite limit, e(oo) <C l/^fi 
we obtain e(L||) = 1/L 2 . Thus, a finite L\\ cuts off the 
low-e tail in the spectral density p(e). 

As an example of the resulting finite-size effects in ther- 
modynamic quantities we now discuss the dependence 
of the susceptibility on Ln. Within the strong-disorder 
renormalization group, it can be found by running the 



renormalization group to the scale ft l = V^f- Beyond 
that scale, e is not renormalized further down. All sur- 
viving layers now have e >• J- 1 and can thus be treated 
as independent. Using (19), the uniform susceptibility of 
a system of size L|| is consequently given as the sum over 
all layers surviving at scale fit, 

X (J,L||)=nn t (J)^ t (J)/2fi L ■ (28) 

At criticality, 5 = 0, this leads to \ ~ 
Ljj [ln(i|[ /a)] 2 * _1/ ^. We emphasize that \ is the sus- 
ceptibility per volume, so L 2 is not simply a geometric 
factor but indicates the divergence of the susceptibility in 
the thermodynamic limit. In the weakly disordered Grif- 
fiths phase, the same calculation gives (up to logarith- 
mic corrections) a non-universal power-law dependence, 
x „ S v+2^a-4>) L 2 - 2 /\ In the W eakly ordered Griffiths 
phase, we need to take into account that long-range order 
is not possible for any finite Lm. Thus all layers surviv- 
ing at scale will again contribute to the susceptibility. 
In contrast to the weakly disordered Griffiths phase, the 
typical moment of a layer is proportional to its thick- 
ness pn L ~ TO ofif2L where mo is the bulk magnetization. 
In the weakly ordered Griffiths phase, we thus obtain 
x ~ 6 v-v4>tI>l 2 + 2 / z . All of our results for the L\\ depen- 
dence of the uniform susceptibility are completely anal- 
ogous to the corresponding temperature dependencies at 
the quantum phase transition in the random transverse- 
field Ising chain in Ref. 6. They are also compatible with 
finite-size scaling using that 1/L 2 scales like e (or, equiv- 
alently, like a magnetic field H). Other observables can 
be worked out in a similar fashion. 

We now turn to the effects of a finite size in perpen- 
dicular direction, i.e., the effects of a finite number of lay- 
ers in our stack. We expect these effects to be particularly 
important experimentally because growing samples con- 
taining a macroscopic number of layers will often be dif- 
ficult. The origin of finite-size effects in L±_ is that finite- 
size samples do not contain rare regions (strongly-coupled 
layers) beyond a certain thickness or, equivalently, they 
do not contain rare regions with e < Q m i n (L±). 

Within the strong-disorder renormalization group, the 
relation between system size and the cutoff energy scale 
Q m in(L± ) can be worked out using the density of surviv- 
ing layers uq. In a typical sample of size L±, the num- 
ber of layers surviving at renormalization group scale f2 
is given by L±uq. The cutoff scale is thus defined by 
L±nn min = 1. At criticality, this implies 

In^j/n^) ~ L\ (29) 

reflecting the activated character of finite-size scaling in 
perpendicular direction. In the two Griffiths phases, we 
obtain 

n mm ~ \5\-" L- ± z . (30) 

As the first example of the resulting finite-size ef- 
fects in thermodynamic quantities, we consider the 
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magnetization-field curve m(h). To do so, we com- 
pare the field-induced renormalization group cutoff f2/j = 
pn h h and the finite-size cutoff ft m in- As long as Clh > 
Qmin, the finite system size L± has only a negligible ef- 
fect on the magnetization. However, L±, cuts off the 
nonlinear low- field tail of m(h) once £lh < Qmin- At 
criticality, this happens for fields below h min given by 
ln(f2/ jhmin) <~ h\. In the weakly disordered Griffiths 
phase, the nonlinear m(h) curve is cut off below h m i n ~ 
s - vz -^(i-4,) L -z_ In the weakly ordered Griffiths phase, 

the calculation is slightly more involved because we first 
need to resolve the relation between Q/, = /iQ h h using 
Mo h ~ m o n nl where m is the bulk magnetization. We 
finally obtain h min ~ \5\<»W-<>(i+z) L^ (1+z) . 

The zero-field susceptibility of a typical sample of per- 
pendicular size L± can be calculated by summing (19) 
over all layers with the spectral density p{e) cut off 
at e = flmin- Alternatively, it can be estimated by 
(dm/dh)h min - At criticality, the susceptibility diverges 
exponentially with system size, x ~ exp(AL^) with A 
a constant. In the weakly disordered Griffiths phase, we 
find a non-universal power law, \ ~ 5 1+vz+2v ^^~^' ] L Z ^ X . 

Finally, we discuss the finite-size behavior of the per- 
pendicular spin-wave stiffness pz\ in the weakly ordered 
Griffiths phase. In Sec. IV C, we showed that pjr vanishes 
in the thermodynamic limit, because the fixed-point dis- 
tribution of J^- is gapless. A finite perpendicular size 
L± establishes a lower bound for J 1 - in a typical sample. 
From (30) we obtain J^ in ~ L^_ z . Thus, the perpendic- 
ular spin-wave stiffness of a typical sample vanishes as 
pjr <~ L 1 ^ 2 with increasing system size (for z > 1). 



VI. CRITICAL DYNAMICS 

It is well known that dynamic critical phenomena show 
stronger rare region effects and Griffiths singularities 
than the corresponding thermodynamic critical phenom- 
ena at classical phase transitions. In particular, rare re- 
gions dominate the long-time dynamics in a conventional 
classical Griffiths phase 20-23 even though they provide 
only small corrections to the thermodynamics. In this 
section, we therefore study the critical dynamics in our 
randomly layered Heiscnberg magnet. 

The classical Heisenberg model does not have any in- 
ternal dynamics, we therefore add a phenomenological 
dynamics to our system. Here, we focus on the sim- 
plest CciSC, el purely relaxational dynamics correspond- 
ing to model A in the classification of Hohcnberg and 
Halperin. 24 Microscopically, this type of dynamics can 
be realized, e.g., via the Glauber 25 or Metropolis 26 algo- 
rithms. Other dynamical algorithms can be studied using 
similar methods (including model J which describes the 
dynamics of real Heisenberg spins). This remains a task 
for the future. 

To characterize the dynamic critical behavior, we cal- 



culate the average autocorrelation function 

C(t) = J d 3 r(0(r, t)<f>(r, 0)) , (31) 

where <p(r, t) is the order parameter at position r and 
time t. In addition, we also determine the dynamic sus- 
ceptibility x( w )- 

Let us begin by considering the dynamics of a single 
layer with moment fi per site and a renormalized local 
distance e from criticality. Because a single layer cannot 
display long-range order, the correlations decay exponen- 
tially in time. The dependence of the correlation (relax- 
ation) time £t on e can be found following the heuristic 
arguments of Bray. 22 He considered a correlation volume 
£ 2 <~ 1/e which he assumed to be in the magnetic state 
with total magnetization Mo ~ /i£ 2 ~ /Lt/e. The relax- 
ation of the magnetization occurs mainly via diffusion of 
the order parameter vector on a sphere of radius \Mq\ 
due to thermal noise (because there are no energy bar- 
riers in the case of Heisenberg symmetry). The noise 
at different points in space and time adds incoherently. 
Thus, according to the central limit theorem, the change 
in magnetization after time t is AM(t) ~ t 1 ^ 2 (iJ,/e) 1 ^ 2 . 
Defining £ t as the time when 5M(£ t ) ~ Mq, we obtain 

6(e) (32) 

At criticality and in the weakly disordered Griffiths 
phase, \i only provides logarithmic corrections to the 
leading 1/e dependence. 

The same result can also be obtained more formally 
from the single-layer Langevin equation 

- -2r (e + 7 2 q 2 )0(q, t) + ^Hq, t) + r?(q, t) 

(33) 

where ft,(q, t) is a time-dependent magnetic field, 7/(q, t) is 
the usual (^-correlated noise and r fixes the overall time 
scale. To find the autocorrelation function of a single 
layer, we solve (33) for h(q, t) = and insert the solution 
into (31). In the asymptotic long-time limit, T et 3> 1, 
we find 

Q(t) ~ exp(-2r et)/(r ei) , (34) 

in agreement with the heuristic estimate (32). Solv- 
ing the Langevin equation in the presence of a field al- 
lows us to calculate the single-layer dynamic susceptibil- 
ity Xz(q,^) = dm(q,co)/dh(q_,uj). For a uniform field, 
q = 0, this results in 

Xi(u) =n 2 /(2e-iLu/T Q ) . (35) 

After having discussed the single-layer dynamics, we 
now turn to the full system. To find the average auto- 
correlation function at time t, we run the strong-disorder 
renormalization group to the scale fl t = 1/t. All layers 
eliminated during this procedure have correlation times 
£ t <C t and do not contribute to the autocorrelation 
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function. Surviving layers have £t 3> t, they thus con- 
tribute proportional to their moment /i per site, giving 
C(t) ~ nn t Mn t - At criticality, this leads to an ultraslow 
logarithmic decay of the autocorrelation function, 

C{t) ~ p^t/to)]*" 1 ^ (36) 

with to a microscopic time scale. In the weakly disor- 
dered Griffiths phase, the same calculation yields, up to 
logarithmic corrections, a non-universal power law 

C(t) ~ <^+"V>(i-</>)f-i/z . (37) 

The same time dependence also follows from averaging 
(34) over the spectral density p(e). The power-law decay 
(37) is much slower than the stretched exponential found 
in conventional classical Griffiths phases. 20-23 Interest- 
ingly, (36) and (37) are reminiscent of the behavior at 
certain classical non- equilibrium phase transitions with 
disorder. 27-29 

The uniform dynamic susceptibility can be computed 
in an analogous manner. At criticality, we find its imag- 
inary part to behave as 

xV^-Hro/u;)] 2 *" 17 *- (38) 

CO 

In the weakly disordered Griffiths phase, we again obtain 
a power law, 

X» ~ tf"+2^(i-*) w V*-i . (39) 

For the local dynamic susceptibility xioc{u), the 
corresponding relations are (l/w)[ln(ro/w)]^ _1// ^ and 
§v +^(i-4>) w i/z-i at cr iti ca lity and in the Griffiths phase, 

respectively. 

At first glance, the above results for C(t) and x"j oc (cj) 
appear to violate the fluctuation dissipation theorem 
which requires x"ioc{ UJ ) = (u>/2T)C(w) where C(ui) is the 
Fourier transform of the autocorrelation function. The 
reason for the disagreement is that the relaxation time 
(32) diverges for the largest rare regions (which corre- 
spond to effective layers with the smallest e). Thus, the 
layers that dominate the long-time tail of C (t) are not in 
equilibrium, and the fluctuation-dissipation theorem is 
not applicable. Technically, the disagreement is caused 
by the fact that (34) cannot be used for the layers with 
the smallest e at any finite time. 

VII. DISCUSSION AND CONCLUSIONS 

In summary, we have investigated the phase transi- 
tion in a three-dimensional randomly layered classical 
Heisenberg magnet. We have employed a strong-disorder 
renormalization group to show that the critical point is 
of unconventional infinite-randomness character. Some- 
what surprisingly, the critical behavior can be found ex- 
actly, making our system one of the very few examples 
of three-dimensional systems with exactly known critical 



exponents. The critical point is accompanied by strong 
power-law Griffiths singularities (which are often called 
quantum Griffiths singularities because they generically 
occur in quantum systems but not in classical systems). 
In addition to the thermodynamics, we have also studied 
the critical dynamics within model A of the Hohenberg- 
Halperin classification. It is characterized by an ultra- 
slow relaxation of the magnetic correlations at criticality 
as well as in the Griffiths phase. 

Our findings can be related to a broader 
classification 9 ' 14 of phase transitions with quenched 
disorder. This classification is based on the effective 
dimensionality of the defects or, equivalently, the rare 
regions. Three classes can be distinguished: (i) If 
the defect dimensionality is below the lower critical 
dimension d~ of the problem, the resulting critical 
point is conventional, and the Griffiths singularities are 
exponentially weak, (ii) If the defect dimensionality is 
exactly equal to the lower critical dimension, the critical 
point is of infinite-randomness type and accompanied 
by power-law "quantum" Griffiths singularities, (iii) 
Finally, if the defects are above the lower critical 
dimension, individual regions can order independently, 
leading to a smeared transition. 

The randomly layered Heisenberg magnet falls into 
class (ii) because the dimensionality of the planar de- 
fects is two, identical to the lower critical dimension of 
the classical Heisenberg model. The results of this pa- 
per are therefore in complete agreement with the general 
classification. It is worth noting, that the behavior of 
a randomly layered Ising magnet is very different. The 
lower critical dimension of a classical Ising model is one, 
thus planar defects are above the lower critical dimen- 
sion. Consequently, the phase transition in a randomly 
layered Ising magnet is smeared by the disorder. 30,31 

We emphasize that the above classification also helps 
resolve the puzzling question posed at the end of Sec. 
IIIB, viz., why systems as different as the randomly lay- 
ered Heisenberg magnet and the McCoy- Wu model (or, 
equivalently, the random transverse-field Ising chain) end 
up in the same universality class. The crucial point is 
that even though these two systems have different order 
parameter symmetries and dimensionalities, the defects 
are exactly at the lower critical dimension in both cases: 
d~ = 2 for the classical Heisenberg model and = 1 for 
the Ising model. These arguments demonstrate why both 
the McCoy- Wu model and our randomly layered Heisen- 
berg model end up having infinite-randomness critical 
points and thus the same scaling scenario; they do not 
yet explain why the two systems share the same critical 
exponent values. The agreement of the exponent values 
follows from the fact that both systems are random in one 
direction which leads to coarse graining in one dimension 
within the strong-disorder renormalization group. More- 
over, the renormalization group fixed point only depends 
on the multiplicative structure the recursion relations (9) 
and (10) and not on model-dependent prefactors. 

Our explicit calculations have been performed in the 
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large- N limit of the order parameter field theory (6). 
However, the critical fixed point stays valid for all TV > 2 
including the case of Heisenberg symmetry, N — 3. To 
see this, we need to confirm that the recursion relations 
(9) and (10) remain unchanged for any N > 2. The 
multiplicative structure of the recursion (9) for J fol- 
lows directly from second order perturbation theory and 
is thus the same for all N. In contrast, the multiplicative 
structure of the recursion relation (10) for the renormal- 
ized distance e from criticality follows from the fact that 
a single layer of a continuous symmetry order parame- 
ter (N > 2) is exactly at the lower critical dimension. 
This implies an exponential dependence of e on the mo- 
ment of the effective layer and thus the multiplicative 
form of (10); for details see the corresponding discussion 
in Ref. 19. Consequently, up to unimportant prefactors, 
both recursion relations remain valid for any N > 2 and 
with them the infinite-randomness critical point scenario 
found in this paper. 

The strong-disorder renormalization group allowed us 
to identify the infinite-randomness fixed point and verify 
its stability. However, it cannot tell whether or not a 
weakly or moderately disordered system will flow toward 
this fixed point. This is due to the fact that for weak 
disorder, the renormalization group recursions are not 
very accurate. To gain further insight, it is useful to look 
at the behavior in the weak disorder limit. The effects 
of weak disorder on a clean critical point are governed 
by the Harris criterion 32 that states that the clean fixed 
point is stable against disorder, if its correlation length 
exponent v fulfills the inequality d r v > 2 where d r is the 
number of dimensions in which there is randomness. In 



our case, d r = 1 and the correlation length exponent of 
the clean 3D Heisenberg model is 33 v « 0.711. There- 
fore, the clean 3D Heisenberg critical point violates the 
inequality d r v > 2 implying that it is unstable against 
weak planar disorder. (It is also unstable against linear 
disorder, d r — 2 (sec Ref. 34), but stable against the 
usual point disorder, d r = 3.) Within a renormalization 
group approach this means that weak planar disorder 
initially increases under renormalization suggesting that 
our fixed point may control the critical behavior for all 
bare disorder strength. A more complete answer to this 
question will likely come from experiment and computer 
simulations. 

Experimental verifications of infinite-randomness 
critical behavior and the accompanying power-law 
"quantum" Griffiths singularities have been hard to 
come by, in particular in higher-dimensional systems. 
Only very recently, promising measurements have 
been reported 35 ' 36 of the quantum phase transitions in 
CePdi-^Rh^ and Nii-^Vj;. We hope that our work 
opens an alternative avenue to observe these phenomena 
in systems that may be easier to study experimentally. 
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